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We report the recent progress in the studies of D-wave heavy quarkonia production from gluon 
fragmentation at the Tevatron and in Z° decays, and the D-wave charmonia production in B 
decays. We show that the color-octet contributions are of crucial importance to the production 
^\ ' of D-wave quarkonia. Signals, as strong as that for J/ip and ip' , should be seen for the D-wave 

charmonium in the processes. These will provide a test of the color-octet production mechanism. 
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Studies of heavy quarkonium production in high energy collisions provide important infor- 
mation on both perturbative and nonperturbative QCD. In recent years, a new framework 
for treating quarkonium production and decays has been advocated by Bodwin^Braaten 
and Lepage in the context of nonrelativistic quantum chromodynamics (NRQCD) £1. In this 
approach, the production process is factorized into short and long distance parts, while 
the latter is associated with the nonperturbative matrix elements of four-fermion operators. 



Qh! This factorization formalism gives rise to a new production mechanism called the color-octet 

mechanism, in which the heavy-quark and antiquark pair is produced at short distance in a 



color-octet configuration and subsequently evolves nonperturbatively into physical quarko- 
nium state. The color-octet term in the gluon fragmentation to Jjipi^P') has been considered 
to explain the J/tp(tp') surplus problems discovered by CDFai. Taking the nonperturba- 
tive (Og // ^( 3 S'i)) and (0% ( 3 5i)) as input parameters, the CDF surplus problem for J/ip 
and tp' can be explained as the contributions of color-octet terms due to gluon fragmenta- 
tion. In the past few years, applications of the NRQCD factorization formalism to J/ip(ip') 
production at various experimental facilities have been studieda. 

Even though the color-octet mechanism has gained some successes in describing the 
production and decays of heavy quark bound systems , it still needs more effort to go 
before finally setting its position and role in heavy quarkonium physics. (For instance, the 
photoproduction data from HERA til put a question about the universality of the color- 
octet matrix elements!, in which the fitted values of the matrix elements (O^^ ^ Sq)) and 
(Og // ^( 3 Pj)) are one order of magnitude smaller than those determined from the Tevatron 
data!). 

In this report, we discuss the production of .D-wave heavy quarkonia under the frame- 
work of NRQCD factorization formalism, including gluon fragmentation at the hadron col- 
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liderll in the Z° decays B, and the D-wave charmonia production in B decays All these 
results show that the color-octet mechanism is crucially important to D-wave charmonium 
production, and the color-octet contributions are found to be over two orders of magnitude 
larger than the color-singlet contributions. This because, in essence, the above processes are 
associated with the color-octet 3 S*i intermediate configuration, which is the same as that in 
J ftp or tp' production at the Tevatron via gluon fragmentation. Therefore, the production 
of D-wave quarkonia would provide a crucial test of the color-octet production mechanism. 
On the other hand, even if the color-singlet model predicts the D-wave production rates too 
small to be visible, they could now be detected after and only after including the color-octet 
production mechanism. It is the color-octet mechanism that could make it possible to search 
for the D-wave heavy quarkonium states, and then complete the studies of charmonium and 
bottomonium families. 



2 Gluon Fragmentation to 3 Dj Heavy Quarkonia at the Tevatron 

We choose a special process to study the gluon fragmentation to color-singlet 3 Dj quarko- 
nium, 

Q*^Qg*; g* ^ 3 Dj gg . (l) 

The decay widths of a virtual quark Q* to color-singlet quarkonium state s Dj by gluon 
fragmentation can be evaluated via 

s 

r(Q* -> Qg*;g* -+ 3 Dj gg) = J dfi 2 T(Q* -> Qg*{fi)) • P{g* -> 3 A/ gg), (2) 



where s is the invariant mass squared of Q*; jx is the virtuality of the gluon, and its minimum 
value squared \J? min = 12mg corresponding to the infrared cutoff as discussed below; P is 
the decay distribution defined as 

P{g* - AX) = — 3 T(g* - AX). (3) 

The calculations of T(g* — > AX) are lengthy but straightforward!^. So, the fragmentation 
functions can be calculated 

n (v , dT(Q* -+*Dj gg Q)/dz 

D 9 ,^ Dj (z,2m Q ,s) = r(Q . - Qg) , (4) 

where z = = 2 — x\ — X2- At high energy limit, the interaction energy s goes up to 
infinity, then the definition of z here is identical with that in Ref.[2] multiplied by a factor 
of two and the fragmentation functions decouple from any specific gluon splitting processes, 
which just reflects the universal spirit of fragmentation. The fragmentation function of 
Eq.dH) is evaluated at the renormalization scale 2mQ, which corresponds to the minimum 
value of the invariant mass of the virtual gluon. 
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The calculation of color-octet fragmentation functions in g* — > 3 Dj( 3 S\,8) processes is 
trivial. They may be obtained directly from color-octet g* — > J/i/j( 3 Si,8) process □, 

D g ^ 3Dj (z,2m Q ) = V °^ q) 5(1 - zKdf'fSi)). (5) 

For the numerical calculations, we choose 

m c = 1.5 GeV, m b = 4.9 GeV, a s (2m c ) = 0.26, a s {2m b ) = 0.19, 

l^(cc)(°)| 2 = 0.015 GeV 7 , |-R(' 6 6)(0)| 2 = 0.637 GeV 7 . (6) 
And then, for color-singlet gluon fragmentation we obtain 
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For color-octet gluon fragmentation, the fragmentation probabilities are proportional 
to the nonperturbative matrix elements (O g Dj ( 3 Si)) which have not been extracted out 
from experimental data, nor from the Lattice QCD calculations. Based upon the NRQCD 
velocity scaling rules, here we tentatively assumeQ 

(ol D2{cg> ( 3 Sx)) « (Of ( 3 ^i)) = 4.6 x 10- 3 GeV 3 (8) 
and further extend this relation to the bb system 

(ol^^fSi)) » (Oj'( 3 ^i)) = 4- 1 x !0 -3 GeV 3 . (9) 

Note that the NRQCD scaling rules tell us that the two matrix elements in (8) (or (9)) are 
of the same order, but not necessarily equal, therefore the assumed relations (g) and @ 
certainly possess uncertainties to some extent. However, from the calculated results below 
we are confident that it will not destroy our major conclusion. From the approximate heavy 
quark spin symmetry relation, we have 

(o^fsi)) « lidpfsi)) « folTsi)) (io) 

for both bb and cc systems. 

Using Eqs.(5), (a), and (S)-(JlQ), we readily have 
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Comparing ( jTT| ) with (j7|), we come to an anticipated conclusion that at the Tevatron the 
gluon fragmentation probabilities through color-octet intermediates to spin-triplet D-wave 
charmonium and bottomonium states are over 2 ~ 4 orders of magnitude larger than that 
of color-singlet processes. As a result, the production rates of 3 Dj charmonium states are 
about the same amount as tp' production rates. Compared with the tp' production at the 
Tevatron, the gluon fragmentation color-octet process plays an even more important role in 
the 3 Dj quarkonium production, and it also gives production probabilities larger than the 
quark fragmentation process t3. 



3 In Z° Decays 

For the D-wave heavy quarkonium production in Z° decays, the dominant process is Z° — > 
3 Djqq. This is similar to Z° — > ipqq discussed in Ref.fll]. Here q represents u, d, s, c or b 
quarks. In the limit m q = 0, we readily have 

TiZ^Djqq) = F(Z _ qq f^ W ~ ?) ~ W 



36 m 
+ [2Li 2 {-^— )-2Li 2 ' ' 



21n(l + £) ln^ + 31n£ + ln 2 ^! + £) 2 }, (12) 
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where Li 2 (x) = — J dt ln(l — t)/t is the Spence function. The calculation with physical 
o 

masses, say, e.g. rrif, = 5 GeV, has also been performed by us, which does not show much 
difference from the case of m q = 0. 

From Eq.(|T2]) (with slight modification due to nonzero m q ), we can get the branching 
ratios of Z° -^ 3 Djqq. Summing over all the quark flavors (q = u, d, s, c, b) with their 
physical masses, we obtain the fraction ratios 

T(Z° ^ 3 D iq q) 4 T(Z° -^ 3 D 2 qq) d T(Z° ^ 3 D 3 qq) d , N 

-\— ^ = 2.0 x lO" 4 , -\— - Q ^ = 3.4 x HT 4 , -\— ^ = 4.8 x 10" 4 . 13 

T(Z U -> qq) T(Z U -» qq) T(Z U -> qq) 

The color-singlet processes include the quark fragmentation contributions and the 
gluon fragmentation contributions. For quark fragmentation, the branching ratios of 3 Dj 
production in color-sing let processes are 2.3 x 10~ 6 , 3.6 x 1CT 6 , 1.7 x 1(T 6 for J = 1, 2, 3, 
respectively. The gluon fragmentation processes are more complicated. For in the most 
important kinematic region the virtual gluon is nearly on its massshell, 3 Dj production 
in the gluon fragmentation color-singlet process may be separated into Z° — ► qqg* and 
g* Djgg. The decay widths of Z° to color-singlet charmonium state 3 Dj by gluon 
fragmentation can be evaluated via 

Aff 

r(Z°^qqg*;g* ^ 3 Djgg)= J dfi 2 T(Z° - qqg*)P(g* - 3 Djgg), (14) 

fi 2 - 
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where the cutoff A = m c is transformed into a lower limit on = 12m^. Summing over 
all the flavors q (q = u,d,s,c,b), we obtain the fraction ratios for the production of 3 Dj 
quarkonium are 4.3 x 1CT 7 , 2.1 x 1CT 6 , 1.2 x 1CT 6 for J = 1, 2, 3, respectively. 



4 D-wave Charmonia Production in B Decays 

We first calculate the color-singlet contribution to 3 D± charmonium production in b decays. 
The matrix elements (0| (cc) v-A 1 3 £*i) can be calculated asi, 

m * )M _ ^ J§. (15) 



So the partial width for color-singlet 8\ production in b decays Q will be 



r(b^cc?DP) + 8,d-+Si + X) = ^ ( ° 8 ^f l)) (2C + -C'_) 2 (l + ^)fo. (16) 
The color-octet contributions to 5j production are similar to that to the J/ ifi production 







<«* + ., + = { ^l + ^£lpi ) (17 ) 



Q1U2 

x(C + + C_) 2 (l + 1? f)f . (18) 



The predicted branching fractions for D-wave charmonium states are 

B(B ^5l + X) = 0.28%, B(B ->6% + X) = 0.46%, 5(5 ^6% + X) = 0.65%. (19) 

We can easily find that the relative production rates predicted above are 51 : 5^ • <5jj — 3 : 
5:7. If we do not take into account the color-octet mechanism the relative rates would be 
SI : 5§ : 5% = 1 : : 0. 

Comparing the predicted production rate of 2 Z)-wave charmonium ( |i"9| ) with that of 
if)', we can see that the former has the same amount as the latter. A similar discussion on 
the D-wave charmonium production in B decays was also presented in Ref.[12]. 



5 Conclusions 



In conclusion, we have shown that in the production of D-wave heavy quarkonia, the color- 
octet processes are always dominant. Their contributions are orders of magnitude larger 
than the color-singlet contributions. According to the NRQCD velocity scaling rules, the 
production rate of D-wave quarkonium is the same order as that of ip' . Among the three 
triplet states of D-wave charmonium, 5% is the most prominent candidate to be discovered 

a Here, we denote the physical D-wave charmonium state as Sj, and 3 Dj represents the leading part of 
the Fock state compenents of 5}. 
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firstly. It is a narrow resonance, and the branching fraction of the decay mode J/^tt + tt is 
estimated to beta 

B{5 C 2 -» JHm + K~) « 0.12, (20) 

which is smaller than that of B{ip' — > J/ipTr + n~) by only a factor of 3. Therefore, 2 D- 
wave charmonium should be observable at the Tevatron via color-octet gluon fragmentation, 
and at the CLEO via B decays. This would provide a crucial test for the color-octet 
production mechanism. We expect the experimental data will be soon obtained. 
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